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ABSTRACT 


A significant disruption of the Quasi-Biennial Oscillation (QBO) occurred 
during the Northern Hemisphere (NH) winter of 2015-16. Since the QBO 
is the major wind variability source in the tropical lower stratosphere and in- 
fluences the rate of ascent of air entering the stratosphere, understanding the 
cause of this singular disruption may provide new insights into the variability 
and sensitivity of the global climate system. Here we examine this disruptive 
event using global reanalysis winds and temperatures from 1980-2016. Re- 
sults reveal record maxima in tropical horizontal momentum fluxes and wave 
forcing of the tropical zonal mean zonal wind over the NH 2015-16 winter. 
The Rossby waves responsible for these record tropical values appear to orig- 
inate in the NH and were focused strongly into the tropics at the 40 hPa level. 
Two additional NH winters, 1987-88 and 2010-11 were also found to have 
large, tropical lower stratosphere, momentum flux divergences; however, the 


QBO westerlies did not change to easterlies in those cases. 
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1. Introduction 


The Quasi-Biennial Oscillation (QBO) consists of downward descending easterly and westerly 
zonal wind regimes that dominate the zonal mean wind variability in the tropical lower strato- 
sphere (100-10 hPa, ~18—30 km in altitude) with a varying (~28 month) period (see Baldwin 
et al. 2001, and references therein). The QBO has been a persistent characteristic of the tropical 
lower stratosphere since observations began in 1953. However, a significant disruption of the QBO 
occurred during the Northern Hemisphere (NH) winter of 2015—16 (Newman et al. 2016; Osprey 
et al. 2016) and several features of this singular disruption imply that a different mechanism may 
have been responsible for the disrupting accelerations than the vertically propagating waves re- 
sponsible for the QBO. Most noticeably, anomalous easterly accelerations occurred in the center 
of the QBO westerlies, a region of weak vertical wind shear, rather than in the strong vertical wind 
shear regions as has been typically observed. 

Vertically propagating equatorial waves are believed to be the principal forcing mechanism of 
the QBO (Lindzen and Holton 1968). Selective filtering of vertically propagating waves by the 
QBO wind distribution coupled with the tendency of the waves to break or thermally dissipate, 
deposit momentum, and thereby dissipate in regions of the QBO wind shear produce appropri- 
ately signed zonal wind accelerations that effectively lower the shear regions by approximately 
1 km month~!. Thus the strength of the wave forcing determines the QBO period. The waves 
responsible are a mix of global scale eastward-propagating Kelvin waves, westward-propagating 
equatorial Rossby-gravity waves and smaller-scale eastward- and westward-propagating gravity 
waves, all originating in the troposphere (Holt et al. 2016). Even relatively small zonal accelera- 


tions can build strong equatorial winds over time as the lack of the Coriolis force at the equator 
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enables the winds to continue in the direction of the acceleration rather than turning as at mid- 
latitudes. 

In contrast to the typical downward propagation of the QBO, based on wave-induced accel- 
erations in the regions of vertical wind shear, Newman et al. (2016) and Osprey et al. (2016) 
found easterlies developing in the region of strong westerlies. Examination of the tropical zonal 
momentum budget by Osprey et al. (2016) showed that the divergence of the horizontal EP flux 
component (Eliassen-Palm flux, see Andrews et al. 1987, page 128) was responsible for the anony- 
mous easterly acceleration near 40 hPa that characterized the 2015—16 disruption of the QBO and, 
in addition, that these EP flux vectors propagated into the tropics from the Northern hemisphere. 
The upward and equatorward EP flux pattern noted by Osprey et al. (2016) is typical of Rossby 
wave propagation in the winter stratosphere (Hamilton 1982), however the effect of Rossby waves 
on the equatorial winds has previously been considered to be small based on idealized model ex- 
periments that showed Rossby waves interacting with the edges of the QBO westerly jet but not 
changing the magnitude of the jet (O’Sullivan 1997). Given the structure of the anomalous QBO 
evolution observed during 2015-16, the potential of Rossby waves to significantly affect the QBO 
needs to re-examined. 

Another possible QBO disruption mechanism would be barotropic instability in the equato- 
rial region. Shuckburgh et al. (2001) showed extensive regions of potential barotropic instability 
associated with QBO westerlies. The relatively small vertical scale of the anomalous easterly ac- 
celeration, centered on ~40 hPa, suggests that barotropic instability may be working to reduce the 
latitudinal wind shear in the region of strong westerlies. In addition to wave forcing we consider 
the possibility of these local wind shear instabilities. 

To characterize the wave forcing responsible for the disruption of the QBO we examine the 


Rossby wave equatorial momentum forcing during the 2015—16 NH winter using global reanalysis 
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winds and temperatures from 1980-2016. This extends the analysis of Osprey et al. (2016) by 
placing the 2015-16 momentum forcing in the context of a 36 year reanalysis climatology. We 
will also examine the possibility of barotropic instability at 40 hPa during the 2015—16 NH winter. 
After describing the data sets used and the analysis procedure (Section 2), we present the mean 
equatorial momentum fluxes and their divergences along with the evolution of the zonal mean 


zonal wind (Section 3), followed by a summary and discussion of the results (Section 4). 


2. Data and Methods 


For this study we use output collections from the Modern-Era Retrospective analysis for Re- 
search and Applications-Version 2, MERRA-2 (Bosilovich et al. 2015) including three-hourly 
instantaneous output on model levels (GMAO 2015b) and monthly averages on constant pres- 
sure levels (GMAO 2015c). The model levels are approximately one kilometer apart in the lower 
stratosphere with ~14 levels between 100 and 10 hPa. In the stratosphere, the pressure levels 
are [100, 70, 50, 40, 30, 20, 10, 7, 5, 4, 3, 2, 1] hPa. MERRA-2 begins in January 1980 and is 
ongoing. The stand-alone MERRA-2 model component generates its own QBO, based on both 
resolved waves and parameterized gravity wave drag (Molod et al. 2015; Holt et al. 2016), thereby 
reducing reliance on observations for the assimilated QBO (Coy et al. 2016). Time altitude cross 
sections of the MERRA-2 QBO zonal mean zonal winds from 1980-2012 are shown in Kawatani 
et al. (2016). Note that all equatorial averages here are based on a 10°S—10°N latitudinal average 
except for Fig. 4 that is based on averages over 5°S—S°N for direct comparison with Osprey et al. 
(2016, their Fig. 2b). 

A QBO composite from MERRA-2 was generated based on the date of the change from zonal 
mean easterlies to westerlies at 30 hPa. The zonal mean zonal winds from the 3 hour collection 


were averaged over a day and from 10°S—10°N before selecting the composite dates of the wind 
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sign change. The composite QBO averages different times of year so that the annual and semi- 
annual cycles tend to average to zero, however, the specific years examined, 2014-16, have both 
annual and semi-annual cycles present. To compare without the annual and semi-annual cycles, 
the monthly averages over the years 1980-2014 were removed when constructing the deviation of 
2014-16 from the composite (Fig. 1c). This procedure mainly removed a semi-annual signal at the 
upper levels shown along with a smaller annual signal. The standard deviation of the composite 
(Fig. 1d) was multiplied by a factor of \/2 to estimate the amplitude of the variability. 

The Eliassen-Palm flux vectors (EP flux, see Andrews et al. 1987, page 128) are a function 
of Rossby wave wind and temperature covariances. The EP flux divergence accelerates the zonal 
mean zonal wind. For this study the EP flux was calculated using the monthly averaged MERRA-2 
data collection. These contain the meridional heat and momentum fluxes (v’T’ and u'v’ where uv’, v’, 
and T’ are zonal wind component, meridional wind component, and temperature respectively and 
the prime denotes a deviation from the zonal mean) needed for the EP flux calculation. However, 
the vertical momentum flux, uw’ (where w is vertical velocity), is not included in the monthly 
averaged collection, so monthly averages of u'w’ were calculated from the 3-hourly assimilation 
output on constant pressure levels (GMAO 2015a). Plotting the EP flux vectors can be problematic 
as they decrease in amplitude at upper levels and in the tropics. To address this issue they are 
plotted only over a limited altitude (70 hPa and above) and latitude (30°S-30°N) range at the 
MERRA-2 constant pressure levels (see above). 

We also used MERRA-2 fields from the monthly mean momentum budget files (GMAO 2015d) 
to distinguish between the parameterized gravity wave drag (GWD) accelerations needed to obtain 
a QBO in the MERRA-2 system (Molod et al. 2015) and the resolved dynamical acceleration, the 
sum of the dynamical and data analysis forcing. These values are accumulated at each time step 


and provide a breakdown of the exact momentum budget. In addition we calculated the monthly 
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averaged zonal mean zonal momentum forcing by the horizontal and vertical EP flux components 
and the residual mean circulation (5°S-5°N) as in Osprey et al. (2016) based on the 3-hourly 
assimilation output on constant pressure levels. 

Also included for February are monthly averaged EP flux vectors and EP flux divergence, nor- 
malized by their standard deviations. As the horizontal component of the EP flux vector is ~2 
orders of magnitude greater than the vertical, a combination of the horizontal and vertical standard 
deviations (horizontal + 100 x vertical) is used to normalize both components, preserving the vec- 
tor directions. The factor of 100 is the order of magnitude of the ratio of the buoyancy frequency 
to the Coriolis parameter at mid-latitudes (N/f,). Since they are normalized by the climatology 
they highlight interannual variability in the flux. 

Along with the EP flux vector, we examine the heat and momentum fluxes separately. Since the 
tropical momentum and heat fluxes are generally an order of magnitude smaller than their winter 
middle latitude values and decrease with altitude, we have normalized these fluxes by their local 
standard deviations when comparing their relative values during individual years. The monthly 
averaged heat and momentum fluxes (GMAO 2015c) were first zonally averaged and then the 
mean and standard deviations were calculated at each latitude and vertical level over the MERRA- 
2 period (1980-2014, 36 or 37 monthly averaged values). After subtracting the multi-year monthly 
mean, the fluxes were then divided by the monthly standard deviation for each location, providing 
normalized values in terms of the local standard deviations. 

The response of the mean meridional circulation to the disrupted QBO was examined by cal- 
culating the residual mean meridional circulation and plotting the vertical component, w*, using 
the same data sets as in the EP flux calculation described above. To focus on the perturbation 
the multi-year monthly average values (Dec 1981 — Feb 2015) were subtracted from each month 


before averaging for the winter season (Dec 2015 — Feb 2016). 
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To assess the possibility of barotropic instability we calculate the meridional gradient of the 


potential vorticity field (Andrews et al. 1987, Eq. 5.3.4): 


oe (acoso )o a (Ppof?_ 


where Q is the Earth’s rotation frequency, a is the Earth’s radius, u is the zonal and time average 
of the MERRA-2 monthly averaged zonal wind component, po is the basic state density, z is the 
log pressure vertical coordinate, and @ is latitude. Note that this differs slightly from the insta- 
bility parameter in Shuckburgh et al. (2001), where only the meridional gradients were examined 
(barotropic instability). Our results showed little contribution from the term involving the vertical 
derivatives (baroclinic instability) so that in this case the barotropic component of the instability 


requirement (q, < 0) dominates. 


3. Results 


The 2015-16 QBO was highly disrupted from its normal behavior. Figure 1 illustrates the time 
height structure of the MERRA-2 zonal mean zonal wind (Fig. la). The longitudinally dependent 
MERRA-2 winds, when zonally averaged, agree well with the local radiosonde winds shown in 
Newman et al. (2016, Fig. 1a) and the zonally averaged assimilation winds presented in Osprey 
et al. (2016, Fig. 1a). The typical zonal wind pattern descent is interrupted by anomalous easterlies 
developing at 40 hPa in early 2016 along with the striking ascent of the westerly winds that began 
in late 2015. In comparison, the composite of the past 14 MERRA-2 QBO cycles (Fig. 1b) shows 
the typical descending shear zones. As in the longer radiosonde record (Newman et al. 2016) the 
MERRA-2 zonally averaged means show that the duration of the QBO westerlies at 40 hPa and 


easterlies at 10 hPa were approximately half of their typical duration. 
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The 2015-16 QBO anomaly with respect to the composite (Fig. lc, the difference between 
Figs. la and b, with the annual and semi-annual cycles removed) shows the vertical extent and 
timing of the QBO disruption. The easterly anomaly at 40 hPa develops over the Nov 2015 —Apr 
2016 period along with the nearly simultaneous development of the westerly anomaly at 10 hPa. 
Note that the rapid appearance of the anomaly at all altitudes (a change over 15 km in altitude 
within a month) is much faster than the usual QBO descent rate (1 km month7!), another indica- 
tion that the 2015-16 dynamics differ from the typical QBO dynamics. The standard deviation of 
the 14 QBO cycle composite (Fig. 1d) shows that most of the QBO variability usually occurs in the 
downward progressing shear zones in agreement with Pawson et al. (1993). Thus the downward 
westerly shear zone in 2014 and early 2015 shows expected variability, while the Dec 2015 and 
later anomaly pattern occurs in regions of weak vertical wind shear and generally low variability, 
indicating an unexpected perturbation of the QBO. 

Figure 2 shows the total zonal mean zonal momentum budget broken down into the parameter- 
ized GWD (red curve) and the resolved dynamics (blue curve). The NH 2015-16 resolved easterly 
accelerations have the largest magnitudes seen during the MERRA-2 period, peaking at -6 m s~! 
month ~! in February 2016. In contrast, the acceleration due to the GWD parameterization, usually 
active during easterly accelerations, peaks at about -2 m s~! month ~! in March and April 2016, 
only about one quarter of its typical value. These parameterized GWD accelerations are positive 
or very small during the months of the anomalous easterly acceleration, November 2016-February 
2016, and contribute little to the momentum budget. This is because the vertical wind shear at 40 
hPa is very small during these months and the parameterization is designed to act strongly in wind 
shear regions. Only after the anomalous easterlies form, creating vertical wind shear near 40 hPa, 


did the GWD parameterization begin to contribute to the zonal momentum budget. 
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Some of the anomalous resolved easterly accelerations were produced by Rossby waves propa- 
gating into the equator from the NH (Osprey et al. 2016). Rossby wave activity propagation from 
the NH into the tropics is proportional to the negative of the horizontal momentum flux (—u’v’, 
see Andrews et al. 1987, chapter 5). Figure 3 shows the time series of the 10°S—10°N, 40 hPa 
monthly averaged horizontal momentum flux (red curve) for the MERRA-2 period. The largest 
peak is seen in the Dec 2015—Feb 2016 period. The Feb 2016 peak is about 50% greater than the 
Jan 2011 maximum. The Dec 2015 and Jan 2016 values are approximately the same as the Jan 
2011 peak. Thus, the NH 2015-16 40 hPa level had the greatest horizontal momentum flux wave 
observed in the 35-year MERRA-2 period. 

As shown by Osprey et al. (2016) the divergence of the horizontal component of the EP flux dur- 
ing November 2015—February 2016 led to the historic easterly acceleration of the QBO westerlies 
at 40 hPa. Fig. 3 shows the monthly averaged 10°S—10°N horizontal momentum flux divergences 
or wind acceleration tendencies (blue curve) during the MERRA-2 period, where negative values 
contribute to a negative EP flux divergence and a negative, or easterly zonal wind acceleration. 
The large amplitude negative peak corresponds to Feb 2016, where there were large momentum 
fluxes (red curve) and an easterly acceleration of the equatorial winds (gray curve). As with the 
momentum fluxes, the Feb 2016 peak is the largest seen at 40 hPa over the 35-year MERRA-2 
period. Comparing with Fig. 2 shows that the horizontal momentum flux divergence is equal to 
about half of the total zonal mean zonal wind acceleration during November 2015—February 2016. 
This implies that the remaining half of the MERRA-2 momentum budget is due to the combination 
of vertical momentum flux divergence and zonal mean circulations since the GWD parameterized 
accelerations are small during the disruption (Fig. 2). 

Different analyses provide an opportunity for comparing their representation of the tropical 


zonal mean momentum budget during the QBO distruption. Here we use a + 5 degree latitudinal 
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average and examine the same momentum budget terms for MERRA-? as presented in Osprey 
et al. (2016, their Fig. 2b) for ECMWE (European Centre for Medium-Range Weather Forecasts). 
Four terms of the 40 hPa, zonal mean momentum budget for Nov 2015 through Apr 2016 are 
plotted in Fig. 4. They consist of the horizontal and vertical EP flux divergence as well as the 
horizontal and vertical residual mean advection. As in Osprey et al. (2016), the horizontal EP flux 
divergence produces the greatest easterly acceleration, peaking in Feb 2016, while the residual 
mean advection terms are relatively small. While the time behavior is similar, the magnitude of 
the Feb peak (~4.5 ms~! month™!) is smaller than in Osprey et al. (2016, ~7.5 ms~! month7'). 
In addition, the MERRA-? vertical EP flux divergence remains small throughout the period shown, 
whereas the Osprey et al. (2016) results show larger values in Mar-Apr 2016. The vertical resolu- 
tion differences between the two analysis system (with ECMWE having higher vertical resolution) 
may contribute to these differences in resolved wave momentum divergence. The missing resolved 
momentum in MERRA-2? is replaced by the GWD parameterization and the analysis increments 
so that the total momentum budget shown in Fig. 2 accurately reflects the changing zonal mean 
zonal wind. 

The NH winter season (Dec-Feb) momentum flux divergence is examined in more detail in 
Fig. 5. The momentum flux divergence tends to be greater during NH winters with QBO wester- 
lies (Fig. 5a). Three winters show exceptionally large magnitudes, 1987-88, 2010-11, and 2015- 
16, with 2015-16 being the greatest. The 1987-88 and 2010-11 NH winters show a weakening 
followed by a strengthening of the QBO westerlies; however mean easterlies do not develop in 
those winters, only during 2015-16. Like the 2015-16 NH winter, 1987-88 coincided with ENSO 
(El Nifio Southern Oscillation), however, the 2010-11 NH winter was about a year after an ENSO. 
Figure 5b further breaks down the winter season into months and shows that, while corresponding 


months in other winters showed some with greater magnitudes, the seasonal average divergence 
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magnitudes were greatest in 2015-16. For comparison, the most recent past westerly QBO NH 
winter, 2013-14, had momentum flux divergence values the were only about one third of the 2015- 
16 magnitudes. 

The mean flow changes can be traced backward to the subtropics using EP flux vectors. This 
wave propagation can be seen in the monthly mean winds and EP fluxes for the 2015-16 winter 
in Fig. 6. In November the equatorial QBO westerlies are centered at about 40 hPa with easterlies 
above. The November EP flux arrows show waves propagating into these westerlies, and across 
the equator — a pattern that is not atypical for QBO westerlies. However, as shown in above, 
the momentum flux divergence is much stronger than in any of the previous westerly phases. 
December shows wave propagation across the equator and the start of small easterly perturbation 
intruding toward the equator. During the Jan—Feb period the westerlies are split into two maxima 
with development of easterlies at 40 hPa with February (Fig. 6d) showing a EP flux pattern similar 
to that found in Osprey et al. (2016). In March the easterlies are fully developed, and continue 
to increase their vertical extent. By April, easterlies completely surround the separated upper 
westerly jet. In summary, during the Nov-Feb period the average lower stratospheric EP fluxes 
extended from north to south across the equator as expected for planetary waves propagating from 
the NH to the SH. A complete understanding of theses waves and their relatively large contribution 
to the momentum budget and flux ( Figs. 4 and 5) needs further investigation. 

Figure 7 illustrates the latitude structure of the horizontal momentum flux, the horizontal mo- 
mentum flux divergence, and the meridional gradient of potential vorticity at 40 hPa for Jan 1998— 
Sep 2016. This figure corresponds to the similar fields shown in Shuckburgh et al. (2001) for the 
30 hPa level. The horizontal momentum flux (Fig. 7a) shows large horizontal momentum flux 
values extending from 30°N across the equatorial region during 2015-16, the time of the anoma- 


lous easterly acceleration. Other years show variability in the strength and equatorial extent of 
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the annual cycle of momentum flux at 30°N with the 3 m*s~? 


contour also extending close to 
the equator during 2010—11 consistent with the large average momentum flux values seen for that 
winter (Fig. 2). The zonal mean zonal wind forcing created by the 2015-16 horizontal compo- 
nent of the momentum flux divergence ( Fig. 7b) shows a corresponding strong region of easterly 
acceleration at the equator extending into the Southern Hemisphere at the time of the anomalous 
easterly acceleration. Note that the 2010-11 westerlies show a northward displacement (but not a 
reversal) of the latitudinal extent of the westerlies during the time of the second greatest equatorial 
horizontal momentum flux values in the MERRA-?2 record (Fig. 2). The potential of the flow for 
instability, gy ( Fig. 7c), shows negative regions typically at the start of the westerly phases but not 
during the anomalous easterly acceleration of 2015-16. Note that the larger wind meridional zonal 
wind shears associated with the beginning of the 2015 QBO westerlies and the time of maximum 
instability are apparent in Newman et al. (2016) their Fig. 2b, a plot of zonal mean zonal wind as 
a function of latitude and time, and furthermore, that these wind shears are greatly reduced at the 
start of the anomalous easterly acceleration. 

Wave activity in the tropics was much higher during the 2015-16 QBO than during the recent 
2013-14 QBO, where the 2013-14 winter provides a more typical example of tropical horizontal 
momentum flux divergence (Fig. 5). The increased wave activity in 2015 compared to 2013 is 
illustrated in Fig. 8, a plot of EPV at 40 hPa averaged over December. The same mean climate 
EPV field has been subtracted from both years to highlight the perturbations. From about 15°S to 
30°N, southwest to northeast sloping, EPV anomalies are seen during 2015 (Fig. 8a) while 2013 
shows smaller amplitude, more zonally oriented EPV anomalies. The zero of the 40 hPa zonal 
mean zonal wind at this time is located at ~15°S so the 2015 EPV orientations are consistent with 
positive momentum fluxes in the region of westerlies. Note that the SH vortex lasted late into Dec 


2015 as denoted by the low EPV anomaly near the South Pole. 
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While all the 2015-16 NH winter months had average or above average tropical momentum 
fluxes, the values for February 2016 were especially notable. Figure 9 shows the local standard 
deviation normalized momentum and heat fluxes at 40 hPa as a function of latitude. The range of 
the previous Februaries (1980-2014) is given by the gray shading. The February 2016 momentum 
flux (Fig. 9a) is nearly 10 standard deviations above the climatology at 10°S. The next largest value 
is in 1983 at nearly 4 standard deviations, much less than the 2016 value. The 2016 momentum 
flux values are greater than 5 standard deviations from 20°S—15°N. As with the momentum fluxes 
the 2016 heat flux (Fig. 9b) stands out from the other years with only 1983 showing an equal 
peak value at 20°N (gray shading). Note that the 2016 heat fluxes are mainly positive north of the 
equator and negative south of the equator indicating upward wave propagation (vertical EP flux 
vectors) in both hemispheres. 

Figure 10 shows February normalized momentum fluxes as a function of latitude and pressure 
for four selected years: 2016 (disrupted QBO), 2014 (a recent more typical westerly QBO), 2011, 
and 1988 (the two years with large amplitude tropical horizontal momentum flux divergence). The 
large tropical values during 2016 are strongly focused at the 40 and 30 hPa levels with values 
greater than 9 standard deviations. February 2016 also shows relatively large positive values (>3) 
at 30°N and 100 hPa. The comparison year, 2014 (Fig. 10b), shows positive fluxes at 40 hPa in 
the tropics; however, they are much smaller (<2) than the 2016 values, and most of the domain 
shows negative values. As in 2013-14, during 2010-11 westerlies continued throughout the winter, 
including February 2011 (Fig. 10c), however, February 2011 resembles 2014 more than 2016 with 
tropical momentum fluxes at 40 hPa peaking near 2 standard deviations. February 1988 (Fig. 10d), 
like 2015-16, was concurrent with a strong ENSO event along with westerlies in the equatorial 
lower stratosphere and the Feb 1988 tropical values are relatively large, peaking at over 2 standard 


deviations, though smaller than the Feb 2016 values. Overall, the 2014, 2011, and 1988 Februaries 
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show negative momentum fluxes at 30°N and 100 hPa, in contrast to 2016. Note that February is 
past the peak month of equatorial horizontal momentum flux divergence for the comparison years 
(Fig. 5). Examination of corresponding plots for December and January (not shown) showed 
horizontal momentum fluxes as large as 3 standard deviations in the lower stratosphere during 
January 2014 and January 2011, and as large as 2 standard deviations in December 1987. These are 
similar to the peak values in found December 2015 and January 2016. None of the corresponding 
positive upper tropospheric values are greater than ~2 standard deviations. Thus February 2016 
especially stands out for its strong horizontal momentum flux values in the NH upper troposphere 
and tropical lower stratosphere. 

Figure 11 compares the February heat fluxes for the same four years. The largest values (-5 
to 4 standard deviations) are found in 2016 at 50 hPa in the tropics. As at 40 hPa (Fig. 9b), the 
field generally switches sign across the equator indicating a strong upward EP flux component 
over most of the tropics. There are also stronger positive and negative values during 2016 in 
the Northern Hemisphere upper troposphere (20-60°N, 150 hPa) than is seen in the other three 
years. Fig. 11 suggests that the tropical waves during 2016 are stronger than average, even in 
the Southern Hemisphere lower stratosphere. While not significant in the MERRA-2 momentum 
budget (Fig. 4), the vertical divergence of EP flux (dependent on the meridional heat flux) in the 
tropics at 40 hPa is shown by Osprey et al. (2016) to be increasing in February 2016 and a leading 
term by March 2016, so that these fluxes may play a role in the later stage of the QBO disruption. 
In addition, the large amplitude meridional heat fluxes seen here in February 2016 suggests that 
the ECMWFE analyses examined in Osprey et al. (2016) can be expected to have correspondingly 
larger amplitude fluxes. 

Figure 12 presents the February anomalous EP flux vectors, again for same four years. Note 


that these are the EP flux vectors normalized by their local standard deviations (Section 2) to 
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highlight the interannual variability and thus differ from the vectors plotted in Fig. 6d. February 
2016 (Fig. 12a) shows larger than average upward fluxes poleward of the Northern Hemisphere 
tropospheric jet (red contours). The large fluxes into the stratosphere turn towards the tropics at 
~40-30 hPa. Large amplitude regions of negative EP flux divergence (red shading) are seen in 
the tropics at those altitudes and in the Southern Hemisphere. In contrast, 2014 ( Fig. 12b) shows 
reduced EP flux into the tropics in the lower stratosphere (poleward arrows). Both 2014 and 2011 
( Fig. 12b and c) show larger than average tropical EP flux vectors, though they are smaller than 
the 2016 case, more upward oriented, and not associated with large anomalous EP flux divergence. 
The 1988 case ( Fig. 12d) shows has anomalous EP flux vectors that are nearly equal in magnitude 
to Feb 2016, however, the tropical divergences are smaller than Feb 2016. None of the three 
additional Februaries examined in Fig. 12 show the large amplitude negative EP flux divergence 
values found in 2016. 

Along with strong tropical wave activity throughout the 2015-16 winter, there was an especially 
large amplitude tropical wave breaking event during early February 2016. The NH polar winter 
of 2015-16 was extremely cold in December and the polar vortex planetary waves were relatively 
weak until late January. The 2015-16 winter then had a very early major final warming event in 
early March (Manney and Lawrence 2016). As the polar planetary wave activity increased in late 
January and a wave breaking event occurred, the tropics responded with an associated strong wave 
event. The exact origin of this strong tropical wave event likely involves some combination of 
stratospheric wave breaking and direct tropospheric forcing that we plan to investigate in future 
modeling studies. Figure 13 shows the evolution of this feature in EPV on the 530 K potential 
temperature surface at 5 day intervals. The winter polar vortex (red shading) displayed a strong 
wavenumber 2 pattern on 31 January 2016 (Fig. 13a) that interacted with the tropical EPV (green 


shading) near 90°E longitude. This produced an intrusion of subtropical air (transparent shading) 
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into the tropics and a wide-in-latitude “knot” of tropical EPV formed and propagated westward 
over equatorial Africa (Fig. 13b). By 10 February (Fig. 13c) the disturbance continued to propa- 
gate westward over the Atlantic Ocean and extended from South American to Africa. While the 
westward propagation slowed somewhat, 15 February found the EPV disturbance centered over 
South America with a long tail of tropical EPV extending south of the equator over the Western 
Pacific. (Note that an animation of Fig. 13, including a comparison with 2013-14, is available as 


supplemental material.) 


4. Summary and Conclusions 


The disruption of the QBO mean zonal wind during the 2015—16 NH winter was associated with 
record strong stratospheric tropical wave activity. This disruption was well captured by MERRA-2 
(Fig. 1). The mean wind disruption was the only event of its kind seen since regular observation 
of the QBO began (Newman et al. 2016). Associated with this record disruption, the tropical wave 
momentum flux at 40 hPa, after very strong values during Dec—Jan, attained a record peak value 
in Feb 2016 (Fig. 3), the largest in magnitude of any month during the 35-year MERRA-?2 period. 
This tropical wave activity was especially focused at the 40 hPa level (Figs. 9 and 10). Initially 
in Nov—Dec 2015, the wave momentum fluxes crossed the equator, reaching the SH easterlies. 
The SH easterlies at 40 hPa then intruded toward and eventually crossed the equator, effectively 
splitting the QBO westerlies (Fig. 6). 


In summary, the boreal winter of 2015-16 showed: 


e record strong momentum and heat fluxes in the tropical lower stratosphere consistent with 


southward and upward wave propagation. 
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e at 40 hPa the developing anomalous easterlies split the QBO westerlies into two distinct 


westerly jets. 


e a large amplitude tropical wave breaking event occurred in February 2016. 


Evidence shown in Osprey et al. (2016) and in Figs. 10 and 12 suggests NH wave generation as 
the most likely source of the anomalous easterly acceleration. However, there is still the question 
of what forced the NH wave generation necessary to cause the 2015-16 QBO disruption. The 
1987-88 and 2010-11 NH winters also showed large tropical momentum flux divergences in the 
tropical lower stratosphere, however, in those years the waves were apparently not of sufficient 
magnitude to reverse the QBO, and westerlies prevailed throughout the winter. So the question re- 
mains about why some NH winters have increased momentum flux divergence and, though some- 
what larger in 2015-16, what specific factors about the 2015-16 winter caused the reversal of the 
zonal mean zonal wind. 

The origins of the 2015-16 NH winter increase in wave forcing needs further investigation. The 
increased wave forcing could have resulted from the naturally large stratospheric-tropospheric 
internal variability, or possibly be tied to specific variability such at that associated with ENSO or 
changed global climate patterns. In particular Newman et al. (2016) (their Fig. 4) showed that the 
tropical upper tropospheric temperatures were much warmer than the MERRA-2 climate record. 
Such warm temperatures may affect tropical and middle latitude wave generation and propagation. 

In the climatological mean, winter season Rossby waves propagate upward and equatorward 
and generally extend into the QBO westerlies. Figure 12 showed that the February 2016 upward 
and equatorward EP fluxes were larger than for the MERRA-2 February average and suggests a 
connection between the middle latitudes and the tropics. However, the heat fluxes for February 


2016 (Fig. 11) showed large values that could be taken to imply more local equatorial Rossby 
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modes as being responsible for the anomalous momentum fluxes, so this possibility is not entirely 
ruled out. However the relatively small contribution of the vertical EP flux divergence to the zonal 
mean equatorial momentum budget (Fig. 4) during the acceleration of the anomalous easterlies 
suggests that the heat fluxes played a relatively small role. We are planning future modeling 
experiments to investigate the specific sources of the anomalous momentum flux. 

Along with the specific cause of the increased wave forcing there remains the need to understand 
why the waves were focused so strongly near 40 hPa in altitude. The QBO westerlies extended 
from ~ 100-5 hPa in the NH fall of 2015, yet the easterly acceleration was strong in a more limited 
vertical region, ~40-30 hPa. This wave focusing allowed the full wave-induced easterly acceler- 
ation to be applied consistently over several months to a relatively confined vertical sub-region 
of the QBO westerlies, adding up to the significant rearrangement of the tropical lower strato- 
spheric winds by the end of March 2016. The intrusion of the easterlies resulting from Rossby 
waves is unexpected given the modeling results of O’ Sullivan (1997) showing only changes in the 
zonal mean wind gradients and not the equatorial jet maximum, so more modeling investigation is 
needed to understand these acceleration. 

Another possibility is a baroclinic, barotropic, or inertial instability associated with the west- 
erly QBO jet. The negative regions of gy of Shuckburgh et al. (2001) suggest the possibility of 
barotropic shear instability associated with the QBO jets. However, the regions of negative q4 
are mainly associated with the increasing QBO westerlies when the meridional wind shears are 
largest. Figure 7 showed that gy was positive during the anomalous easterly acceleration making 
instability of the large scale flow unlikely in this case. Moreover, the mean instability would need 
to be maintained over the several months that characterized the anomalous easterly acceleration. 

More detailed diagnostic and model forecast studies are needed to resolve meridional circulation 


changes associated with this 2015-16 disrupted QBO and to test the ability of seasonal forecast 
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systems to encompass and predict such a disruption of the QBO. As noted by Newman et al. 
(2016) and Osprey et al. (2016) the normally downward propagating westerlies showed an upward 
propagation (or displacement) in 2016 at altitudes above ~30 hPa in the lower stratosphere (Fig. 1). 
Figure 14 plots the Dec 2015—Feb 2016 vertical component of the residual mean circulation (with 
multi-year means removed), w*. The calculated w* field shows upward motion above ~40 hPa 
centered at ~5°S. The upward values of ~1 km month™! are the same order of magnitude as the 
observed upward displacement and suggest that the meridional circulation response to the easterly 
acceleration at 40 hPa played a role in the observed upward displacement. The upward progression 
of the westerlies can therefore be expected to modify the transport and distribution of stratospheric 
trace gases and aerosols. 

The 2015-16 disruption of the QBO provides an opportunity for improving forecasting in the 
tropical lower stratosphere, especially on seasonal time scales, as it provides a specific example 
of how the QBO responds to changes in wave forcing. In this context the winters of 1987-88 
and 2010-11 provide additional examples of strong wave momentum forcing that lacked the zonal 
wind reversals, so that any forecasting improvements should encompass these winters as well. 
Along with developing the ability to forecast a major disruption of the QBO, the QBO disruption 
of 2015-16 may require re-evaluation of the normally high QBO seasonal prediction skill (Scaife 


et al. 2014). 
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2, zonal mean zonal wind acceleration due to parameterized gravity wave drag (red, m 
s-!month~!) and the resolved dynamics (blue, m s-!month~!), and zonal mean zonal wind 
(gray, x10 ms7~!). Vertical lines denote the start of a year; b) expanded time coordinate 
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Also shown is the zonal mean zonal wind (m s~!, right y-axis, black curve). Note that the 
labeled acceleration units of 3, 6,9 m s~! month7! correspond to 0.1, 0.2, 0.3 m s-! day 
Monthly mean momentum flux divergence (m s~!month~!) for NH winter (December, Jan- 
uary, and February) plotted as a) a function of the 10°S-10°N zonal mean zonal wind and 
b) as a function of year. The two digit labels denote the January year. In panel b the winter 
average (wide bars) is broken down into the three monthly averages (narrow bars) where 
blue denotes easterly and red denotes westerly zonal mean zonal winds. 


Monthly averaged zonal mean zonal wind plotted from 30°S-30°N and from 200-4 hPa for 
the months: a) Nov 2015, b) Dec 2015, c) Jan 2016, d) Feb 2016, e) Mar 2016, and f) Apr 
2016. Westerlies are yellow-red and easterlies are green-blue with 5 m s~! contours. Also 
plotted are the EP flux vectors (blue arrows) at 70 hPa and above. . 


Latitude time contour plots at 40 hPa of a) the horizontal momentum flux (m*s~7), b) the 
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Fig. 13. 


Fig. 14. 


February zonally averaged momentum flux for a) 2016, b) 2014, c) 2011, and d) 1988 as 
function of latitude and pressure. The values are non-dimensional in terms of standard devi- 
ations over the years 1980-2014 with a contour interval of one standard deviation. Negative 
values are shaded gray. The red horizontal line denotes the 40 hPa level. 


Same as Fig. 10 for heat flux. 


February zonally averaged zonal wind (10 ms~!, red contours, positive values gray shaded) 
for a) 2016, b) 2014, c) 2011, and d) 1988 as function of latitude and pressure. The arrows 
denote normalized EP Flux deviations from the 1980-2014 February climatology. They are 
normalized as described in Section 2 and plotted so that 5 degrees of latitude corresponds 
to 1 standard deviation. The red (blue) filled regions denote negative (positive) EP Flux di- 
vergence anomalies (non-dimensional, standard deviations, 0.5 contour interval, white con- 
tours). The filled contours start at +1.5. The blue horizontal line denotes the 40 hPa level. 
38 


EPV on the 530 K potential temperature surface for 00 UTC on a) January 31, b) February 5, 
c) February 10, and d) February 15 of 2016. The green colors denote values from ~ -15-15 
PVU, red denote values >100 PVU, and purple denote values <-50 PVU. Latitude lines at 
-60, -30, 0, 30, and 60 degrees. Longitude lines at -135, -90, -45, 0, 45, 90, and 135 aie 
The 530 K surface is approximately at 40 hPa near the equator. ; 


The vertical component of the residual mean circulation (km month~!) averaged Dec 2015 


— Feb 2016 as a function of latitude and pressure. The multi year (Dec 1980- Feb 2015) 


monthly means have been subtracted. Negative values are shaded. 
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FIG. 6. Monthly averaged zonal mean zonal wind plotted from 30°S-30°N and from 200-4 hPa for the months: 


a) Nov 2015, b) Dec 2015, c) Jan 2016, d) Feb 2016, e) Mar 2016, and f) Apr 2016. Westerlies are yellow-red 
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and above. 
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FIG. 7. Latitude time contour plots at 40 hPa of a) the horizontal momentum flux (m?s~7), b) the divergence of 
the horizontal momentum flux (ms~!day~'), and c) the meridional gradient of potential vorticity (10~'!'m~!s~!). 
The black contours highlight the +3, -0.1, and 0 contours in a,b, and c respectively. The green curves denote the 
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606 Fic. 9. Zonally averaged momentum (a) and heat (b) fluxes at 40 hPa for February 2016 (red curve, 10°S- 
so7_1O°N) and plotted as functions of latitude. The values are non-dimensional in terms of standard deviations over 


se the years 1980-2014. The gray shaded regions denotes the February normalized range over 1980-2014. 
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Fic. 10. February zonally averaged momentum flux for a) 2016, b) 2014, c) 2011, and d) 1988 as function of 


latitude and pressure. The values are non-dimensional in terms of standard deviations over the years 1980-2014 


with a contour interval of one standard deviation. Negative values are shaded gray. The red horizontal line 


denotes the 40 hPa level. 
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Fic. 11. Same as Fig. 10 for heat flux. 
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Fic. 12. February zonally averaged zonal wind (10 ms~!, red contours, positive values gray shaded) for 
a) 2016, b) 2014, c) 2011, and d) 1988 as function of latitude and pressure. The arrows denote normalized 
EP Flux deviations from the 1980-2014 February climatology. They are normalized as described in Section 2 
and plotted so that 5 degrees of latitude corresponds to | standard deviation. The red (blue) filled regions denote 


negative (positive) EP Flux divergence anomalies (non-dimensional, standard deviations, 0.5 contour interval, 


white contours). The filled contours start at +1.5. The blue horizontal line denotes the 40 hPa level. 
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Fic. 13. EPV on the 530 K potential temperature surface for 00 UTC on a) January 31, b) February 5, c) 


February 10, and d) February 15 of 2016. The green colors denote values from ~ -15-—15 PVU, red denote 
values >100 PVU, and purple denote values <-50 PVU. Latitude lines at -60, -30, 0, 30, and 60 degrees. 
Longitude lines at -135, -90, -45, 0, 45, 90, and 135 degrees. The 530 K surface is approximately at 40 hPa near 


the equator. 
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624 FIG. 14. The vertical component of the residual mean circulation (km month~!) averaged Dec 2015 — Feb 


es 2016 as a function of latitude and pressure. The multi year (Dec 1980— Feb 2015) monthly means have been 


ss subtracted. Negative values are shaded. 
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